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The universal form for the average scattering intensity from systems undergoing order-disorder 
transitions is found by numerical integration of the Langevin dynamics. The result is nearly identical 
' for simulations involving two different forms of the local contribution to the free energy, supporting 

the idea that the Model A dynamical universality class includes a wide range of local free-energy 
forms. An absolute comparison with no adjustable parameters is made to the forms predicted by the 
theories of Ohta-Jasnow-Kawasaki and Mazenko. The numerical results are well described by the 
former theory, except in the cross-over region between scattering dominated by domain geometry 
and scattering determined by Porod's law. 
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I. INTRODUCTION 



e 

*^ I Phase ordering by quenching from a region of the phase diagram where a material is uniform to one where several 
<**] ■ phases coexist at equilibrium provides an important technique for creating multi-phase materials that are inhomo- 
geneous on mesoscopic length scales. Since the macroscopic properties of such materials can be quite different from 
those of the constituent phases and depend sensitively on the mesoscopic structure, a solid understanding of pat- 
tern formation in phase-ordering systems is important to several branches of materials science. Examples include 
t-H , precipitation strengthening in metals and fabrication of glasses Q . 

^ ■ Universality allows one to describe the structure and dynamical properties of diverse phase-ordering systems using 
models that only take into account properties such as conservation laws and order-parameter symmetries. Two 
important universality classes involve only local relaxational dynamics and an order parameter that can be represented 
' as a scalar field The order parameter is not conserved in the first class, called Model A. This can be used, for 

example, to model anisotropic magnets and alloys undergoing order-disorder transitions. In the second class, called 
Model B, the order parameter is a locally conserved quantity, and relaxation proceeds by diffusion away from regions 
ON . of high chemical potential. When hydrodynamic modes and strain effects can be ignored, binary mixtures and alloys 
■ are described by this model. 

Scaling is the hypothesis that the behavior of the system over a large range of length scales can be described 
H i in terms of a single characteristic length R. A necessary condition for scaling is that R must be well separated 
from any other microscopic or macroscopic length scales present in the system. For many phase-ordering processes, 
the characteristic length has been found to have a power-law dependence on the time, r, elapsed since the quench, 
R ~ t" . In non-conserved systems the power law is readily observable with n = 1/2. For conserved systems the 
q . late-time growth exponent is n = 1/3, although processes occurring at early times may mask this behavior. Both of 
these universality classes are important. However, the work reported here is restricted to non-conserved systems, for 
which more detailed theoretical results are available. 

Here, scaling and universality in phase ordering are tested by comparing numerical solutions of the Langevin 
equation involving two different forms for the local part of the free-energy functional: the ip 4 Ginzburg-Landau form 
and a piecewise-linear triangular form. The numerical details are presented in Sec. II. The agreement between the 
structure factors found for both forms is discussed in Sec. III. Combined with a mapping to sharpen the interfaces, 
numerical solution of the Langevin equation allows comparison to analytic theories for the structure factor, without 
fitting parameters. In Sec. IV three theories are reviewed and compared to the numerical results. The Ohta-Jasnow- 
Kawasaki theory agrees with the simulations at small and large wavevectors with noticeable deviations only for 
intermediate wavevectors. The theories by Mazenko are qualitatively correct, but the zeroth-order theory agrees with 
the simulations better than the second-order theory. Sec. V is a brief summary of our results. 
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II. NUMERICAL MODEL 



A general model of phase ordering can be constructed from a free energy composed of a local term with two 
degenerate minima and a nonlocal term representing the contribution from spatial fluctuations || , 

Fmr,r)}= J dr|/[^(r,T)] + i|V^(r,r)| 2 | , (1) 

where the scalar order-parameter field t/>(r, t) depends on position r and time r. The dynamics of the model are 
governed by a Langevin equation, 

—d^ = - W(r,r) + ^ T ^ 

where the first term on the right-hand side corresponds to deterministic relaxation towards a minimum of the free 
energy, and the second term represents thermal fluctuations. The thermal noise is assumed to be Gaussian with zero 
mean and correlations given by 

(r,(T,r)r,(v'y))=S(r-T')S(r-T'). (3) 

The strength of the noise is given by the normalized temperature e, which is the only parameter in the model 
after rescaling the order parameter, space, and time |^|. The primary effect of thermal fluctuations is to introduce 
randomness at early times ||. In the late-time scaling regime the dynamics of this model are controlled by a zero- 
temperature fixed point, and thermal fluctuations can be ignored. 

The local part of the free-energy functional T is usually chosen to have the Ginzburg-Landau form, 

/[^(r,r)]=-^ 2 (r,r) + ^ 4 (r,r), (4) 

which is a double well with degenerate minima. The corresponding dynamical equation is 

^Elll = (i + v 2 ) ^(r, r) - ^ 3 (r, r) + ^(r, r) , (5) 

which is commonly called the time-dependent Ginzburg-Landau (TDGL) equation. 

To test universality for non-conserved scalar order-parameter systems, comparisons need to be made with a model 
different from the standard TDGL model. The new model still needs to follow the non-conserved dynamics of Eq. (||) 
for a scalar order parameter with two degenerate equilibrium values. One obvious aspect to change is the form of 
the local part of the free-energy functional, f[ip(r, r)]. A choice quite different from the Ginzburg-Landau form is a 
piecewise-linear triangular double well, which is not harmonic near the minima. This form is compared to the ip 4 
Ginzburg-Landau potential in Fig. 1. It has been chosen to match the Ginzburg-Landau potential at the extrema, 
i.e. /(±1) = —1/4 and /(0) = 0. The Langevin equation for this potential is 

f -1/4, V<-1 
tyjr, t) r ,,^2,, v J +1/4, -1<^<0 

— = ^(r,r)+V#,r)- _^ < <A < +1 " ^ 

I +1/4, ip>+l 

The simulations with both potentials were conducted on square lattices with periodic boundary conditions and 
lattice constant Ar = 1 . The system size considered was L x = L y = L = 1024. The Laplacian was approximated 
by an 8-point form [Q, and and a simple Euler integration scheme with Ar = 0.05 was used to collect data every 25 
time units up to a maximum of r = 2000. Results for the TDGL equation were averaged over 100 realizations; only 
10 realizations of the piecewise-linear model were needed to confirm its agreement with the TDGL results. 

III. UNIVERSALITY 

The average structure of the system at a specific time after the quench can be described by the order-parameter 
correlation function 
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C(|r-r / |,r) = (^(r,r)V(r , ) r)) 



(7) 



where (• • •) denotes averaging over the ensemble of quenches as well as the volume. When scaling holds, the correlation 
function can be expressed in terms of a time-independent master curve, 

C(r, T )=C(f). (8) 

where f = r/R(r). While microscopy techniques can be used to measure C, scattering experiments probe its Fourier 
transform, S(k, t), called the structure factor. The structure factor is also related to a time- independent form, 

S(k 1 T) = R d ( T )F(q) (9) 

where q — kR{r) and d is the dimension of the system. 

When infinitely sharp domain interfaces are randomly oriented throughout the system, the small-f form of the 
correlation function is linear, C(f) — 1 — af + • • •, where a is proportional to the surface area to volume ratio for the 
domains ||. This corresponds to a power-law decay of the scaled structure factor at large q, F(q) = (2n/a)q~( d+1 \ 
called Porod's law These results are valid for 1/q much smaller than the characteristic length, but much larger 

than the domain interface width. 

The assumption of infinitely sharp interfaces is used in most situations, e.g. experimental and analytic studies. By 
design, the domain interface in the simulations has a width of about We compensate for this by using a non-linear 
mapping Q of ip(r,r) to ±1 before finding the Fourier transform ^(k, r) Q. The structure factor is 

S(k,T)H|^(k,T)| 2 }, (io) 

where ^(k, r)| 2 is the scattering intensity associated with a particular domain pattern ip(r,r), and (• • •) represents 
averaging over the ensemble of quenches. The simulation results used to estimate the structure factor were found by 
"onion-shell" binning of the two-dimensional scattering intensity into a one-dimensional function fl of fc, averaging 
over all trials, and then scaling the result using the average wavevector of each bin. 

Including noise in the simulation has three drawbacks. The first is purely practical: generating Gaussian random 
numbers for r\ is computationally expensive. The second drawback has a physical basis. Ohta |Hj] has shown that 
thermal fluctuations in the TDGL equation retard growth of the characteristic length and cause a broadening of F(q). 
Thirdly, the thermal noise causes the structure factor to cross over to the equilibrium q~ 2 behavior for large q, thus 
masking the Porod tail. For these reasons, we have simulated the TDGL equation at zero temperature with 

the early-time fluctuations implemented by an initial condition where ij}{r, r = 0) consists of independent random 
numbers distributed uniformly on [—0.1, +0.1]. 

On the other hand, thermal fluctuations are essential for simulating the piecewise-linear model. For this model at 
T = 0, the domain growth stops after a finite time. The time until the evolution stops, as well as the corresponding 
characteristic length, depends on the random initial condition. If a system which has stopped coarsening is briefly 
heated, the growth continues for a while and then stops again. In addition, making the potential shallower prolongs 
growth and yields larger domains. The arrested growth appears to result from trapping into metastable configurations 
associated with finite wavelength corrugations on the interfaces. A constant, moderate level of thermal noise com- 
pletely prevents this trapping phenomenon, so noise with e = 0.1 was included in simulations of the piecewise-linear 
model. 

An estimate of F (g) was constructed for each local form by combining the estimates at different simulation times 
with the restrictions r > 500 and k < 0.5. For k larger than this, lattice effects become important. The estimates 
of F(q) for both models are presented in Fig. 2, with R(t) = y/2r (see Sec. IV). The agreement seen in Fig. 2(a), a 
log-log plot of F vs. q, is very good. The Porod tail at large q can be highlighted by plotting q d+1 F(q) against q, as 
shown in Fig. 2(b). This was found independently of the structure factor, using q = kR(r) before the binning and 
averaging. The agreement is quite good well into the power-law tail. The claim that the details of the local part of 
the free-energy functional do not influence the universal behavior appears well supported by the numerical models. 

IV. ANALYTIC THEORIES 

Several theories for the universal form of the scaling function for non-conserved systems exist. They can be compared 
to the numerical integration of the TDGL equation without adjustable parameters. 

One theory was developed by Ohta, Jasnow, and Kawasaki (OJK) ||TQ , |ll|| . The OJK theory starts from a diffusion 
equation for a Gaussian auxiliary field, u(r, r), 
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du(r, ■ 
dr 



= DV 2 u(y,t), 



(11) 



where the interfaces in the inhomogeneous material are defined to be the set of r such that u(r, r) = 0, and the 
order-parameter field is obtained by the mapping ip(r,r) = sign [u(r, r)]. The coefficient D = 4pd is the diffusion 
constant for the interface. The factor pd = (d — l)/d results from assuming that the interfaces are randomly oriented. 
As time progresses, the characteristic size of the single-phase domains grows as -Rojk(t) = V Dt, in accordance with 
the Lifshitz- Allen-Cahn theory of domain growth driven by surface tension [ p~2|JT^ ] . The two-point correlation function 
for the order-parameter field in this model has the simple analytic form 



CoiK(r, t) 



-arcsm 



exp 



(12) 



Since the system is isotropic, the Fourier transform of this correlation function can be written in terms of a radial 
integral over a Bessel function, and the scaling form of the structure factor is flllH 



-FojkO?) 



2 (27r) d / 2 



duw d [exp {w 2 ) - l]- 1/2 (qw) 1 - d/2 J d/2 (qw) 



(13) 



Except for the dimensional dependence of pd in R(r), this result agrees with the perturbation calculation by Kawasaki, 
Yalabik, and Gunton |1 



Another theoretical approach to phase ordering has been developed by Mazenko |16|-fl9|. This theory also starts 
with a Gaussian auxiliary field u(r, r). In contrast to the OJK theory, u(r, t) is interpreted as the distance from 
r to the closest interface, and it is mapped onto the order-parameter field using the equilibrium interface profile 
ijj(r, t) = ip eq (u(r, t)). The characteristic length is then defined in terms of the RMS distance to the nearest interface, 
Ro{t) = y/n(u(r, t) 2 ). Here the subscript is used to denote the zeroth-order theory developed in Refs. [|6| [L8|. 
Using the Gaussian properties of u(r, r), the TDGL equation can be used to find a closed, nonlinear partial differential 
equation which describes the scaling form of the correlation function 



d 2 C (^) 
dv 2 



1 



P.QV 



dC (v) 
dv 



tan 



A(v) 



o, 



(14) 



where v = t/Rq(t) is the scaled length. The characteristic length is i?o( T ) = V^^o T with po an undetermined 
parameter. The correlation function must obey the boundary condition Cq(0) = 1. As a consequence, the third term 
in Eq. jl^ ) acquires a 1/v singularity at v=0. This singularity must cancel against the (d— l)/v factor in the second 
term, resulting in the requirement that dCo(v)/dv = yj2/Tr{d — 1) at v = 0. Therefore, the first derivative of Cq(v) 
at v—Q is not available to be used as the second, independent constant in the solution of Eq. (|l4|). Instead, one must 
determine the value of po such that Cq(v) becomes integrable as v — > oo. For d = 1 the asymptotic form agrees 
with that found analytically for the kinetic Ising model with Glauber dynamics ]20| , pT[ |, while in the limit d — > oo 
the function Co becomes identical with the OJK result flq|. For finite d > 1, determining Co is equivalent to the 
numerical solution of a nonlinear eigenvalue problem ||l7,|22| , [23"|| . 

Recently, Bray and coworkers PJ24| and Yeung et al. [ 25T, have criticized the use of Gaussian fields in constructing 



theories of phase ordering. Noting this criticism, Mazenko |19| has expanded his theory to include a u(r, r) whose 
distribution is an expansion around a Gaussian. The theory with the second-order correction, here marked with the 
subscript 2, produces a new differential equation for the scaled correlation function, C2(i>), 



d 2 C 2 {v) 
dv 2 



P2V 



d-l\ dC 2 (v) 



tan 



H 2 (v) 



{§ [C 2 (v) + H 2 (v)}} 







(15) 



Here H 2 (v) is a new function governed by 



dH 2 
dv 



where C 2 = dC 2 /dv and 



bC' 2 ~ yJ(C 2 ) 2 +a Q q 2 b(l-b) 



b=~H 2 (v)t a n{^[C 2 (v) + H 2 (v)]} 



(16) 



(17) 
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A corrected eigenvalue [i 2 7^ A*o results, and the characteristic length becomes R%(r") = \I^[i 2 t with the scaled length 
v = r/R 2 {r). A new eigenvalue q 2 is introduced, and the resulting double eigenvalue problem must solved numerically 

While experimental tests ultimately determine the validity of any theory, finding the characteristic length a priori 
is difficult for experiments, as well as many simulation techniques. The characteristic length may be taken as a fitting 
parameter, with which data can be fit well to all three theories. Such fitted comparisons for Monte Carlo simulations 
made by O JK jllj showed quantitative disagreements in the tail of the correlation function. Later Oono and Puri |27[] 
used cell dynamics simulations to show that the deviations might be caused by the nonzero interface width present 



in simulations. Blundell, Bray, and Sattler 24 also used cell dynamics simulations to test the theories. To compare 
their simulation results to theory without fitting, they reduced R(t) to a parametric variable by plotting values for 
the correlation function for the order parameter against values for the correlation function for the square of the order 
parameter for many (r, r) . The simulation results plotted this way show good scaling, but the two Gaussian theories 
give the same result for this particular scaling function [p4| 25|. As a consequence, this method does not provide a 
strong test of the two forms, which are, in fact, noticeably different, as we will see below. The question of how well 
these theories predict the scaling function C(f) therefore remains open. 

Numerical simulation of the TDGL equation can be compared to the theories without fitting parameters since 
analytic forms for R(t) are known. We have chosen to work in terms of -Rojk since Rq and R2 depend on numerically 
determined eigenvalues. The TDGL estimates for the two-dimensional F(q) for several times are presented in Fig. 3, 
along with the predictions of the theories. Fig. 3(a) is a log-log plot, showing good collapse of the simulation data. 
Concerning the difference between the theories of OJK and Kawasaki et al. 113], much better agreement with the 
simulations is found using the factor predicted by OJK, so we only consider that theory. In Fig. 3(a) the obvious 
differences between the theories shown are the magnitude at small q and the behavior near the shoulder at intermediate 
q. For both of these features, the simulation data agree significantly better with the OJK theory than with those of 
Mazenko. 

Fig. 3(b) highlights the power-law tail at large q. The scaling of the simulation data is quite good through the large 
peak for all of the times presented here. Where the simulation data scales, it agrees very well with the OJK theory. 
There is a noticeable second peak in the OJK and second-order Mazenko theories that is absent in the zeroth-order 
Mazenko theory. Unfortunately, the simulation data cannot be used to test the presence of a second peak because 
it does not scale in this region for the practically attainable simulation times. Instead, a pronounced trough is seen, 
which decays with time. Simulations for significantly larger systems and longer times would be needed to achieve 
asymptotic scaling in this region. 

At values of q higher than this, the scaling regime also has not been reached, preventing an estimate of the Porod 
amplitude (2n/a) for the simulations. However, a clear trend in the data is apparent, indicating an amplitude which 
is significantly higher than those predicted by both of Mazenko's theories. In fact, the second-order theory has the 
smaller Porod amplitude. This is a result of the relationship |l9| ] 



y/q 2 (q 2 + 2) 

a 2 =a — , (18) 

92 + 1 

where q 2 is the second-order eigenvalue. Since a 2 > ao for any value of q 2 , the predicted Porod amplitude is necessarily 
smaller for the second-order theory. The disagreement with simulations increasing with the order of the theory is also 
true for the decay of the order-parameter autocorrelation function, although estimates for both of Mazenko's theories 
are better than that for OJK 

Very recently, Emmott has presented a perturbation expansion about the OJK form for the correlation function, 
with a lowest-order term of order 1/d 2 P9|| . Our numerical results presented here suggest that the OJK theory 
may provide a good starting point for a more complete theory of the dynamics of phase separation in systems with 
nonconserved order parameter. However, we note that the lowest-order correction provided by Emmott's approach 
would produce a change in the Porod amplitude (which is well described by OJK), while it shares with OJK the 
inability to correctly describe the decay of the order-parameter autocorrelation function [f7 19 2§f|. 



V. CONCLUSIONS 



Numerical simulations of the TDGL equation have been used to confirm the universality of Model A systems with 
different forms for the local free-energy contribution. With no adjustable parameters, the scaling form of the structure 
factor from simulations is found to agree quantitatively with the OJK theory at scaled wavevectors q < 3. At large 
q, the numerical results have not yet entered the asymptotic scaling regime, even at our latest simulation times. 
However, the trend is toward OJK and away from both the Mazenko theories. Although the OJK theory has been 
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found here to describe well the scaling form of the structure factor, it does not correctly describe the decay of the 
order-parameter autocorrelation function fj\jffi$j . 

Mazenko's theories describe the estimated scaling form only qualitatively. Interestingly the zeroth-order theory 
gives better agreement with simulations than does the second-order theory. Since, aside from d — 1, d — 2 is expected 
to be the dimension where the theoretical results are the most different, computational cost has lead us to forgo three 
dimensional simulations at this time. Another aspect of universality not considered here is the effect of anisotropy 
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FIG. 1. Comparison of the two local parts of the free-energy functional, f[tp(r, t)], used in this study: the Ginzburg-Landau 
form (solid curve) and a piecewise linear function (dashed lines). Both potentials give the same estimate for the structure 
factor, indicating that the Model A universality class does not depend on the specific form of /. 
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FIG. 2. Simulation estimates for both the Ginzburg-Landau (O) and piecewise- linear ({>) models of (a) the universal form 
of the structure factor and (b) scaling plot emphasizing the Porod's law tail at large q. The estimates incorporate all data for 
wavevectors k < 0.5 and time r > 500. The scaling forms for the two models agree with each other, supporting the idea that 
dynamical universality is not affected by the detailed shape of the local part of the free-energy density f[ip]. The OJK (solid 
curve), Mazenko zeroth-order (dashed curve), and Mazenko second-order (dotted curve) theories are included for reference. 
The simulations show better agreement with the OJK theory than with the other theories. 
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FIG. 3. Scaling of the structure factor for the two-dimensional time- dependent Ginzburg-Landau model. The scaling variable 
is q — kRojK(T~) . The simulation data for different times scale quite well, using the characteristic length -Rojk(t) = y/Apdj, 
except at large wavevectors where lattice effects are important. The curves represent the theoretical scaling forms discussed in 
the text: OJK (solid), Mazenko zeroth-order (dashed), and Mazenko second-order (dotted), (a) For small values of q (q < 3), 
the simulations agree quantitatively with the OJK theory, (b) Plot emphasizing the Porod's law tail at large q, the trend 
with time indicates that the OJK theory gives a better estimate of this amplitude than either of the Mazenko theories. The 
theoretical amplitude of the second-order theory is lower than that for the first-order theory. 
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